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ABSTRACT 

We present results for neutrino-nucleon cross sections for energies up to 10 21 eV, of 
relevance to the detection of ultrahigh energy galactic and extragalactic neutrinos. 
At the highest energies, our results are about 2.4 times larger than previous estimates. 
Using these new cross sections, we predict neutrino telescope event rates for the upward 
moving muons initiated by the neutrino interactions in the Earth and for the contained- 
vertex events in the PeV range due to neutrino-electron interactions. We show that 
future neutrino detectors, such as AMANDA, BAIKAL, DUMAND and NESTOR have 
a very good chance of detecting neutrinos which originate in the Active Galactic Nuclei. 

INTRODUCTION 

The Active Galactic Nuclei (AGN), with typical luminosities in the range 10 42 to 
10 48 erg/s, are believed to be the most powerful individual sources of radiation in the 
Universe. These extragalactic point sources are also considered as prodigious particle 
accelerators presumably powered by the gravitational energy of matter spiraling in to a 
supermassive black hole, though the mechanism responsible for the conversion of grav- 
itational energy to luminous energy is not presently understood. Recent detection of 
energetic photons (£ 7 ~ 100 MeV) from about 40 AGNs by the EGRET collaboration!!! 
and of TeV photons from Mkn 421, Mkn 5010 and most recently from 1ES2344+514 
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by the Whipple collaboration" have created new excitement in the field of high-energy 
gamma-ray physics. If the observed photons are decay products of 7r°s produced in 
hadronic interactions in the disk surroundingthe AGN, then AGNs are also powerful 
sources of ultrahigh-energy (UHE) neutrinos. □ Unlike photons, which are absorbed by 
a few hundred gm/cm 2 of material, TeV neutrinos have interaction lengths on the order 
of 250 kt/cm 2 and thus can provide a direct window to the most energetic processes in 
the universe. 

The advantage of the long interaction length translates to a challenge in the detec- 
tion of neutrinos. Interaction rates increase with energy, but the fluxes of UHE neutri- 
nos are steeply falling functions of neutrino energy. Cerenkov detection of muons from 
interactions of muon neutrinos in the rock or ice surrounding the detector is feasible.! 
More difficult is the detection of charged-current interactions of electron neutrinos. 
Large-area air shower arrays or large volume underground detectors may be adequate 
for the detection of electron neutrinos, especially near the IF-boson resonance in v e e 
collisions. Theoretical calculations of the neutrino-nucleon and neutrino-electron cross 
sections are instrumental in evaluating event rates for neutrino telescopes. 

Here we present results! for charged current and neutral current cross sections 
for energies up to 10 21 eV obtained using new parton distributions measured in ep 
collisions at HERA.! We also discuss how detection of UHE neutrinos depends on 
these cross sections and on the neutrino fluxes from UHE neutrino sources. Event rates 
for muon neutrino conversions to muons are compared with earlier results based on 
older parton distribution functions.! We also present results for contained events with 
higher threshold energies. 

SOURCES OF UHE NEUTRINOS 

A variety of sources may contribute to the neutrino flux at the surface of the 
Earth. Three types of sources are discussed here: atmospheric neutrinos from cosmic- 
ray interactions in the atmosphere, neutrinos from active galactic nuclei, and cosmic 
neutrinos from extragalactic cosmic ray interactions with the microwave background 
radiation. Model predictions for neutrino fluxes from these three types of sources are 
shown in Figure 1. Atmospheric neutrinos! (ATM), while interesting in their own right, 
mask extraterrestrial sources for E u < 1 TeV. Consequently, we restrict our discussion 
to neutrino energies above 1 TeV. 

The TeV photons observed by Whipple collaboration! may be byproducts of 
hadronic cascades initiated by the protons generated within the AGN accretion disk of 
gas, or in the jets, which interact with matter or radiation in the AGN disk, to produce 
pions whose decay products include both photons and neutrinos. The structure of the 
corresponding hadronic cascade is: 

PP — > 7T + X 
PJ — ► 7T + X 

np — > n + X 
7T° ^7 + 7 

V — > + v e + e 

If charged and neutral pions are produced in equal proportions and photons originate 
in hadronic cascades, simple counting leads to equal fluxes of photons and + v^. The 
flux of v e + v e equals half of the flux of v^ + v^. The observed photon energy spectrum 



Figure 1: Muon neutrino plus antineutrino fluxes at the Earth's surface: angle-averaged 
flux from atmospheric neutrinos (ATM), diffuse flux from active galactic nuclei (AGN-NMB, 
AGN-SP and AGN-SS) and cosmic neutrinos (CR-2 and CR-4). The Frejus upper limit 1 on 
a neutrino flux in excess of atmospheric neutrino flux is indicated at 2.6 TeV. The dotted line 
indicates the vertical flux of atmospheric \x + p, from Ref. 15. 



is a power-law witl 
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for 100 MeV< E 1 < 2 TeV, and the same for neutrinos. We have chosen three 
representative fluxes of neutrinos from AGN, each corresponding to the diffuse flux 
integrated over all AGNs. These fluxes are shown in Figure 1. The Nellen, Mannheim 
and Biermann flux0 (AGN-NMB), which comes from assuming that pp collisions are 
the dominant neutrino source, is parameterized by: 



= 1.13 x 10- 12 (^/TeV)- 2 cm- 2 s- 1 sr- 1 GeV- 

dE v 



with the v e + v e spectrum assumed to be 1/2 of + v^. The neutrino luminosity of 
a source is normalized to the observed diffuse x-rays and 7-rays. The NMB param- 
eterization is valid for E v < 4 x 10 4 GeV. In our calculations described in the next 
section, we have used this parameterization up to E v = 10 8 GeV. A somewhat differ- 
ent assumption of the luminosity is used by Szabo and ProtheroeEH (AGN-SP) in their 
extended model of neutrino sources, yielding a higher normalization of dN/dE v at 1 
TeV. Above E u > 10 6 GeV, the AGN-SP follows a steeper power law, 



dN/dE v ~ E 



-3.5 



which accounts for the lack of protons at even higher energies required to produce 
neutrinos. The Stecker and Salamon flux@ (AGN-SS) contains contributions from 
both pp and pj interactions in the accretion disk and has a nearly constant value of 
dN/dE v up to E v ~ 10 5 GeV. 



Two models of neutrino fluxes from cosmic ray interactions with the microwave 
background" are labeled CR-2 and CR-4 in Figure 1. The fluxes depend on the redshift 
of the cosmic ray sources. Maximum redshifts contributing are z max = 2 and z max = 4, 
respectively. 

The electron neutrino plus antineutrino fluxes, to a good approximation, are equal 
to half of the fluxes shown in Figure 1. 

UHE MUON NEUTRINOS 

The primary means of detection of muon neutrinos and antineutrinos is by charged- 
current conversion into muons and antimuons. The long range of the muon means that 
the effective volume of an underground detector can be significantly larger than the 
instrumented volume. For example, a 10 TeV muon produced by a charged- current 
interaction in rock will propagate several kilometers in water-equivalent distance units 
before its energy is degraded to 1 TeV. 

Backgrounds to AGN sources of z/ M + include atmospheric neutrinos and at- 
mospheric muons. Muons produced by cosmic ray interactions in the atmosphere mask 
astrophysical signals unless detectors are very deep underground, muon energy thresh- 
olds are set very high, or one observes upward-going muons. We evaluate here event 
rates for upward-going muons produced in the rock surrounding the detector, for muon 
energy thresholds above 1 TeV and 10 TeV. 

The neutrino-nucleon cross section comes into the calculation of the event rate 
in two ways. The probability of conversion — > fi is proportional to the vN charged 
current cross section. In addition, the neutrino flux is attenuated by passage through 
the Earth. In the next section we describe our calculation of the neutrino-isoscalar 
nucleon {vN) cross section. The vN charged-current reaction is the dominant source of 
neutrino interactions except in a very narrow energy window at the M^-boson resonance. 

SMALL-i PARTON DISTRIBUTION FUNCTIONS AND a{vN) 
The inclusive cross section for + N — > /i~ + X is given by 

<fa - ' 2G ^ lE " M ™ .[ xq (x,Q>) + x(l- y mx,Q% (1) 
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where x = Q 2 /2Mu, y = u/E v , with — Q 2 the momentum transfer between the neutrino 
and muon, and v the lepton energy loss in the lab frame, v — E v — E^. M is the 
mass of the nucleon and Myy is the mass of the IV-boson, while the Fermi constant 
is Gp = 1.16 x 10~ 5 GeV~ 2 . Taking the target as isoscalar nucleons, in terms of the 
parton distribution functions for the proton, 

q(x, Q 2 ) = + + s s + b s (2) 

q(x,Q 2 ) = ^±^ + c s + t s (3) 

where we have written explicitly valence (v) and sea (s) distributions. 

The general form of the cross section shows that at low energies, where the four- 
Fermi approximation is valid, a ~ E. At higher energies, the W-boson propagator 
plays an important role. The value of (Q 2 ) saturates at ~ Myy, and x ~ M^/(2ME u y) 
decreases. For neutrino energies above 10 5 GeV, the small-x (x < 3 x 10~ 2 ) behavior 
of the parton distribution functions becomes important for the evaluation of the cross 
section. 



Figure 2: Comparison of the light-quark sea at Q 2 = for various parton distribution 
functions. Of the MRS distributions, D_ (A') is the most (least) singular. 



Figure 3: The charged-current cross section for the CTEQ-DIS, CTEQ-DLA, EHLQ-DLA, 
MRS A', MRS G and MRS D_ parton distribution functions. The data point, an average of 
ZEUS and HI, is from Ref. 17. 



Neutrino charged- current interactions have been measured directly in laboratory 
experiments for neutrino energies up to E v = 300 GevJll Charged- current ep scattering 
at HERA, equivalent to E v = 47.4 TeV, can be translated to a value of cr(z/iV).0 Recent 
ZEUS and HI measurements at HERA0 of F 2 ep at small-x (10" 4 < x < 10~ 2 ) and for 
a large range of Q 2 , 4 GeV 2 < Q 2 < 1600 GeV 2 have provided valuable information 
about parton densities at small-x and \ow-Q 2 . To evaluate the neutrino- nucleon cross 
section at ultrahigh energies, extrapolations beyond the measured regime in x and Q 2 
are required. 

There are two main theoretical approaches in the evolution in Q 2 of parton densi- 
ties: Gribov-Lipatov-Altarelli-Parisi@ (GLAP) evolution and Balitskii-Fadin-Kuraev- 
LipatovB (BFKL) evolution. In the GLAP approach, parton distribution functions are 
extracted at modest values of Q 2 and evolved to higher scales. The BFKL approach 
involves a leading a s ln(l/a;) resummation of soft gluon emissions, which generates a 
singular behavior in x at an initial scale Qo, 



xq s (x,Q 2 )~x-°- 5 (4) 

for small x, which persists at higher values of Q. In our extrapolation of the parton 
distribution functions outside the measured region, we use GLAP evolution with input 
at Qo = 1.6 GeV, 

xq s (x,Ql) ~ x~ x . (5) 

The value of A is determined by fits to deep- inelastic scattering and hadron-hadron data 
by the MRS§ and CTEcjl Collaborations. The MRS set A' has A = 0.17, the MRS 
set G has A = 0.07 while the MRS set D_ has A = 0.5. All of the MRS distribution 
function are fitted using the MS factorization scheme. The CTEQ-DIS, using the deep- 
inelastic scattering factorization scheme, has A = 0.33. These distribution functions 
are extrapolated using the power law fit to the distribution functions at x = 10~ 5 
and Q = Mw- We have also extrapolated the leading-order CTEQ distributions using 
the double-log approximation.il For reference, the Eichten et a/M parton distribution 
functions, extrapolated using the double-log approximation, are also shown. The spread 
in values for the parton distribution functions is an indication of the uncertainty in 
evaluating the vN cross section. 

For each of these sets of distribution functions, we have evaluated the neutrino- 
nucleon cross section. Figure 3 illustrates the range of predictions as a function of 
neutrino energy. Also shown is the average of HI and ZEUS effective neutrino nucleon 
cross sect ions. lUl There is excellent agreement among the predictions of the MRS D_, 
G, and A' distributions and the CTEQ3 distributions up to E v « 10 7 GeV. Above that 
energy, our DLA modification of the CTEQ3 distributions gives a lower cross section 
than the full CTEQ3 distributions (CTEQ-DIS), as expected from its less singular 
behavior as x —>■ 0. At the highest energy displayed, the most singular (MRS D_) 
distribution predicts a significantly higher cross section than the others. Above about 
10 6 GeV, the EHLQ-DLA distributions yield noticeably smaller cross sections than 
the modern distributions. Plots similar to Figure 3 for antineutrino-nucleon charged 
current interactions, as well as neutral current interactions, can be found in Ref. 6. 
For charged current and neutral current interactions, for 10 15 eV < E v < 10 21 eV, the 
cross sections follow a simple power law, for example 

/ P x 0.402 
36_2 / 



vcc(vN) = 2.69 x 10 _Jb cm 



1 GeV 



Table 1: Number of upward fi + fi per year per steradian for A = 0.1 km 2 and 
E mm = x TeV _ 



Fluxes 


EHLQ-DLA 


CTEQ-DIS 


AGN-SS0 


82 


92 


AGN-NMB0 


100 


111 


AGN-SP0 


2660 


2960 


atmI 


126 


141 



Table 2: As in Table 1, but for E™ in = 10 TeV. 



Fluxes 


EHLQ-DLA 


CTEQ-DIS 


AGN-SS0 


46 


51 


AGN-NMB0 


31 


34 


AGN-SpO 


760 


843 


atmI 


3 


3 



NEUTRINO TELESCOPE EVENT RATE 

In order to calculate the number of upward-moving muons that can be detected 
with neutrino detectors such as AMANDA, BAIKAL, DUMAND II and NESTOR,! 
we fold in the neutrino flux and its attenuation in the Earth with the probability that 
a neutrino passing on a detector trajectory creates a muon in the rock that traverses 
the detector. 

The attenuation of neutrinos in the Earth is described by a shadow factor S(E V ), 
equivalent to the effective solid angle for upward muons, normalized to 2tt: 



dS{E v ) = J_ 
dtt 2vr 



exp 



(-z{9)N a <j vN {E v )), (6) 



where N A = 6.022 x 10 23 mol" 1 = 6.022 x 10 23 cm" 3 (water equivalent) is Avogadro's 
number, and z(8) is the column depth of the earth, in water-equivalent units, which 
depends on zenith angle.il The probability that the neutrino with energy E v converts 
to a muon is proportional to the cross section and depends on the threshold energy for 
the muon E™ m : 

P,{E V , Ef) = a cc {E v )N A {R{E v , £f)>, (7) 

where the average muon range in rock is (R) @ A more detailed discussion appears in 
Ref. 6. 

The diffuse flux of AGN neutrinos, summed over all AGN sources, is isotropic, so 
the event rate is 



Rate = A J dE u P^(E u , E™ n )S(E, 



dN v 



dE, 



V 



given a neutrino spectrum dN u / dE u and detector area A. As the cross section increases, 
increases, but the effective solid angle decreases. 

Event rates for upward muons and antimuons for a detector with A = 0.1 km 2 
for Ef n = 1 TeV and Ef n = 10 TeV are shown in Tables 1 and 2. The CTEQ-DIS 
distribution functions are taken as representative of the modern parton distribution 



Table 3: Downward resonance i> e e — > W~ events per year per steradian for a detector 
with effective volume V e R = 1 km 3 together with the potential downward (upward) 
background from and interactions above 3 PeV. 



Mode 


AGN-SSB 


AGN-SPEI 


W -> v^n 


6 


3 


W — > hadrons 


41 


19 




33 (7) 


19 (4) 




13 (3) 


7(1) 



function sets, and compared with the EHLQ-DLA event rate predictions. The muon 
range is that of Ref. 26. 

The theoretical predictions for ultrahigh-energy neutrinos from AGNs yield event 
rates comparable to, or in excess of, the background rate of atmospheric neutrinos for 
E mm _ 1 Te y_ The AGN _gp rate ig large comparec i to the AGN-NMB rate because 

additional mechanisms are included. Flux limits from the Frejus experiment are in- 
consistent with the SP flux for 1 TeV< E u < 10 TeV.0 The atmospheric neutrino 
background is greatly reduced by requiring a 10 TeV muon threshold, though AGN 
induced event rates are reduced as well. The flatter neutrino spectra have larger con- 
tributions to the event rate for muon energies away from the threshold muon energy 
than the steep atmospheric flux. 

We have evaluated the event rates using the other parton distribution functions 
shown in Figure 2. Event rate predictions are unchanged with the other modern parton 
distributions because all these distributions are in agreement in the energy range E v ~ 
1 — 100 TeV. However, our results for event rates are about 15% larger than for the 
EHLQ structure functions. This is due to the fact that EHLQ parton distributions were 
based on the CERN-Dortmund-Heidelberg-Saclay measurements of neutrino-nucleon 
structure function all which had low normalization of about 15%. 

UHE ELECTRON NEUTRINOS 

Finally we consider event rates from electron neutrino and antineutrino interac- 
tions. For v e N (and v e N) interactions, the cross sections are identical to the muon 
neutrino (antineutrino) nucleon cross sections. Because of the rapid energy loss or 
annihilation of electrons and positrons, it is generally true that only contained-vertex 
events can be observed. Since electron neutrino fluxes are small, an extremely large 
effective volume is needed to get measurable event rates. There is one exceptional case: 
resonant formation of W~ in v e e interactions at E v = 6.3 PeV. The resonant cross 
section is larger than the vN cross section at any energy up to 10 21 eV. In Fig. 4 we 
present neutrino-electron cross sections. 

We note that, at the resonance energy, upward- moving electron antineutrinos do 
not survive passage through the Earth. However, the contained events have better 
prospects for detection. The contained event rate for resonant W production is 

10 r dN- 

Rate = ^V eS N A J dE Pe a Vee {E, e )S{E v - e )^. (9) 

We show event rates for resonant H^-boson production in Table 3. The background is 
for events with E v > 3 PeV. 

From Table 3 we note that a 1 km 3 detector with energy threshold in the PeV 
range would be suitable for detecting resonant z/ e e — > W events. However, the u^N 



Figure 4: Cross sections for neutrino interactions on electron targets. At low energies, from 
largest to smallest cross section, the processes are (i) v e e — > hadrons, (ii) u^e — > fiu e , (iii) 
v e e -> v e e, (iv) v e e —* (v) v e e ->• v e e, (vi) v^e -> u^e, (vii) v^e -> v^e. 

background may be difficult to overcome. By placing the detector a few km under- 
ground, one can reduce atmospheric-muon background, which is 5 events per year per 
steradian at the surface of the Earth for > 3 PeV. 

SUMMARY 

In summary, we find that detectors such as DUMAND II, AMANDA, BAIKAL 
and NESTOR have a very good chance of being able to test different models for neutrino 
production in the AGNsEa For E™ m = 1 TeV, we find that the range of theoretical 
fluxes leads to event rates of 900-29,600 upward-moving muons/yr/km 2 /sr originating 
from the diffuse AGN neutrinos, with the atmospheric background of 1400 events/yr 
/km 2 /sr. For E™ m = 10 TeV, signal to background ratio becomes even better, with 
signals being on the order of 500-8,400 events/yr/km 2 /sr, a factor ~20-300 higher than 
the background rate. For neutrino energies above 3 PeV there is significant contribution 
to the muon rate due to the v e interaction with electrons, due to the VT-resonance 
contribution. We find that acoustic detectors with 3 PeV threshold and with effective 
volume of 0.2 km 3 , such as DUMAND, would detect 48 hadronic cascades per year 
from W — > hadrons, 7 events from W — > fiu^ and 36 events from and interactions 
with virtually no background from ATM neutrinos. 
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